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Abstract 

A new formalism for the optimal control of quantum mechanical physical 
observables is presented. This approach is based on an analogous classical 
control technique reported previously. |E| Quantum Lagrange multiplier func- 
tions are used to preserve a chosen subset of the observable dynamics of in- 
terest. As a result, a corresponding small set of Lagrange multipliers needs to 
be calculated and they are only a function of time. This is a considerable sim- 
plification over traditional quantum optimal control theory. ||] The success of 
the new approach is based on taking advantage of the multiplicity of solutions 
to virtually any problem of quantum control to meet a physical objective. A 
family of such simplified formulations is introduced and numerically tested. 
Results are presented for these algorithms and compared with previous re- 
ported work on a model problem for selective unimolecular reaction induced 
by an external optical electric field. 
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I. INTRODUCTION 



Different approaches and paradigms for controlling molecular motion have been pro- 



posed. [IU| Various results indicate that the final state distribution can be controlled in 
many instances. Manipulating the interference between two or more routes to the same 
(degenerate) final state has been suggested to achieve control of the final state distribution. 

Another approach employs the laser field to guide the wave packet motion utilizing two 
electronic potentials. |7]|],[25[] In addition to these particular schemes, a general optimal de- 
sign formalism for the quantum control problem has been developed. [TT|-|T3[ The capability 
of designing laser pulse shapes with this optimal design formalism has been demonstrated 
for many applications including selective excitations, |p7j^2| selective bond breaking for tri- 



atomic molecules, || control of curve-crossing reactions p9],|30] an< ^ control of the electric 
susceptibility of a molecular gas. |HJ] Furthermore, it is known that multiple control solutions 
will likely exist for any particular system. |33] The multiplicity of solutions gives a range 
of flexibility for field design which is especially important for adaptive feedback laboratory 
control. p],g 

The existence of these multiple optimal solutions provides the freedom to develop sim- 
plified methods to find a control field. This paper builds on this observation to present a 
different approach to quantum optimal control theory. The approach only requires that the 
feedback Lagrange variables be scalars to preserve desired observable expectation values. In 
contrast, the previous generic method @j2^j required the propagation of a typically non- 
linear Schrodinger type equation for the Lagrange functions causing considerable numerical 
complexity. The simplicity of the new approach presented here should allow for the study 
of molecular control of larger dynamical systems. A family of simplified approaches will be 
introduced. The new approaches will be applied to the selective optical dissociation of a 
model triatomic molecule. 

The paper is organized as follows. In section II we present the general quantum control 
dynamics equations and the proposed simplifications. In section III we apply the theory to 
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a selective unimolecular reaction. The computational method to solve the control dynamics 
equations is presented in section IV. In section V numerical results for these methods are 
presented and compared to previous results. We discuss and summarize the results in section 
VI. 

II. THE QUANTUM CONTROL PROBLEM AND ITS SIMPLIFICATION 

In this section we present the quantum control dynamics equations for a molecule where 
the control is a laser field. The optimal control theory seeks a field pulse to steer the 
molecular motion from the original state to achieve a desired final objective. Consider the 
Hamiltonian H to have the form 

H = H mol (z) + H mt (R,e(t)) (1) 

where H mo i(z) is the Hamiltonian of the undisturbed molecule, H int (FL, e(t)) represents the 
field-molecule interaction, and e(t) is the electric field vector, z is defined as an operator 
vector z = [z±, zi---Zin\ = [R, P] = [Ri, Rn, Pi, Pn] containing the coordinate R and 
congugate momentum P operators of the system. N is the number of degrees of freedom in 
the molecule. In the dipole model, the interaction Hamiltonian is 

H mt (R,e(t)) = -fi(R)-e(t) (2) 

where /x(R) is the dipole moment vector. The control design formalism seeks the field e(t). 
The molecular motion evolution at time t is described by \ip(t)) which obeys 

,nMSl = „ m)) (3 ) 

with the initial state being 1^(0)). The knowledge of the wave function permits evaluation 
of any dynamical observables, and in particular (z(t)) = (ip(t)\z\ip(t)) . We will assume that 
the control objective can be expressed in terms of (z(t)) although other expectation values 
(e.g., bond energy, etc.) could just as well be treated. 
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A. Conventional formulation of optimal control theory 



In this subsection we will briefly summarize the conventional approach to optimally 
designing fields for molecular control. We will build on this background and some inherent 
design freedom to introduce alternate and simplified formulations in subsection II. B. 

To design the control field e(t) that best achieves the desired objective we introduce the 
cost functional J[(z),e] 

J[(z),e] = $[<z(T)>] + f dt[Z x {(z{t))) +£ 2 (e(t))\ (4) 

J o 

The first part $[(z(T))] is a functional that measures the deviation from the desired physical 
objectives at time T. The second part involves the constraint function £i((z(t))) that takes 
into account any restrictions on the variables (z(t)) over the interval < t < T (e.g., to 
avoid undesired regions of phase space, or products, or in order to help guide the molecular 
system to the desired objective. f22f ) The last part, with the field cost function £ 2 (e(t)), 
expresses the desire to minimize the energy fluence or possibly introduce other biases in the 
field. 

The constraint that the Schrodinger equation be satisfied is assured through introduction 
of the Lagrange multiplier vector |A(t)) along with its complex conjugate. Thus, the full 
cost functional is given by 

J[(z},e} = J[(z} 1 e}-J o T dt(\m^^^-H\m)} 

- £ dt[zh^^ + m)\H}\Ht)) (5) 

The minimization of J[(z), e] leads to a link between the objective functional and the optimal 
solution e(t) supplied by the Lagrange multiplier vector. 

The quantum variational problem is given by 5J[(z),e]=0 which minimizes the cost 
functional J[(z),e] with respect to e(t), \4>(t)) and |A(t)). The variation with respect to 
\\(t)) gives rise to the Schrodinger equation Eq. (|3|). The variation with respect to \i/j(t)) 
leads to the following equation: 
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dl JAii =jr iA( t ) > _g_JLm„ Wt)> ( 6 ) 

This equation for the Lagrange multiplier vector has the same form as the Schrodinger 
equation along with an inhomogeneous term. The final time condition is 

The gradient of the cost functional with respect to the field is 

The solutions of the set of equations Eqs. (|3|) and (§|)-(BD produces the optimal field. 

In the above approach, the unconstrained cost functional J[(z),e] was introduced to 
conserve the Schrodinger equation, and the control equation of motion for the Lagrange 
multiplier vector gives similar dynamics to that of Schrddinger's equation. In this case the 
Lagrange multiplier vector takes into account the state of the molecule at each instant of 
time in the quantum control process. A key observation is that there are generally multiple 
solutions, and possibly a denumerably infinite number, to the control design equations. |33| 
We will take advantage of this flexibility below. The Lagrange multiplier plays the role of 
guiding the dynamics to a particular solution. 

B. Simplified formulation optimal control theory 

The physical cost in equation (|J) only depends on the expectation values and external 
interaction field and not explicitly on the wave function. We can take advantage of this 
observation to simplify the process of achieving control solutions. The time dependence of 
the expectation values is governed by the equations 

which can collectively be written as 



= (iP(t)\fi(z,<t))\m) * = h-,2N (10) 

where the functions fi(z, e(t)) may be readily identified as momentum or coordinate deriva- 
tives of the Hamiltonian. We may now write a new unconstrained cost functional that 
preserves the dynamical equation fllOl), 



(11) 



J[<z),e]= J[(z),e 

T 2N 
i=l 



o — ; dt 

In this unconstrained functional there is an implicit dependence on the Schrodinger equation, 
and this point will become important below. It is significant to note that the functional in 
( p|) retains exactly the same physical objective through J[(z), e] as in equation fl5|). Here we 
only alter the feedback bias towards a particular control solution by the choice of Lagrange 
functions. 

The minimization of J[(z),e] is considered with respect to e(t), \i(t), and (zi(t)). The 
variation with respect to \i(t) gives the quantum equation of motion for the expectation 
values, Eq. ([II]). The variation of \tp(t)) is not explicitly treated, but we do need to consider 
the interpretation of the variation 5(tp(t)\fi(z, e(t))\ip(t)). The various alternate optimal 
control approaches are based on different interpretations for this variation. We introduce 
the variation of {x/j(i)\fi(z,e(t))\i(}(i)), as 

^^(t)|/ t (z,e(t))|^(t))^g^(t)| ^ ( ^ e(t)) ^#(t)) 

j=l OZ 3 

«£(#)l fl/ ' ( f' c(t)) l#)){^(*)) (12) 

where (8zj(t)) is 

(5z 3 {t)) = m)\^m)) (is) 

and \ipit)) satisfies equation (^|). Thus, in employing the functional in Eq. (|TT| ) an im- 
portant point is that the Schrodinger equation for \ij)(t)) is not approximated so that the 
true molecular dynamics is fully retained. Secondly, the original cost functional J[(z),e] 
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in equation (|j) is retained, implying that any control solution e(i) obtained through this 
alternate formulation is just as valid as obtained the conventional way in section II. A. Thus, 
the variation in equation fll^D should be thought of as guiding the design process from one 
valid solution to another equally valid one. The only question at this point is whether this 
new approach can guide the process to achieve designs e(t) that produce quality control. 
The computations in section III will show that the approach can achieve excellent results. 
Considering the above arguments we have the full variation of J[(z), e], as 



T 2V 



,/u i=l 



3=1 

d£i((z(t))) 



(Szi(t)) (14b) 



+ jf *l^ + gWW*)l^hKfl>WO. (-) 

The boundary conditions at time T, 

0$((z(T))) . 

A - (T)= d{ Zl (T)) Z= '-' ' ( } 

are obtained through the Eq. ( |14a| ) . As, we require 5J[(z),e] = 0, the second equation 
derived from Eq. ( |14b| ) is 



T 2N 



j=i i 



and the gradient with respect to the field is 



<M*)> = °> ( 16 ) 



The Lagrange multipliers are chosen to obey the following equations 
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The control dynamics equations to solve the quantum control problem are given by Eqs. 
(§), (|18D and (|17]) with the final conditions for the Lagrange multipliers given by Eq. (|15|). 
The results based on this approach will be called method I. 

The fundamental distinction between solving these equations and those of the standard 
treatment in Eqs. (0), ©-© lies in the nature of the Lagrange multipliers. Equation (|6|) is a 
partial differential equation, while Eq. (]TJ) is a much simplier ordinary differential equation. 
A key point is that the control e(i) generated from the new formulation should be fully 
acceptable physically, despite the simplified form in Eq. (fIB|), since the Schrodinger equation 
is not approximated. The new formulation takes advantage of the multiplicity of solutions 
to steer the design process from one possible field to another. In some cases the control 
might even be similar to that found with the full conventional control formulation (this is 
the case in the examples of Section III). In other cases the control field may be different, 
but importantly retention of the original cost functional assure a proper test of the achieved 
results. 

It is worthwhile to explore if further simplified formulations for the Lagrange multipliers 
can also be successfully introduced. One may replace Eq. (|T8|) by the following: 



y Am ^-((z(t)),e(t)) ^((z(t))) 

~dt h j »<*(*)> ^wr 2N (19) 

This is a classical like formulation for Xi(t), but quantum mechanics is still fully retained 
in solving for \ip{t)) and (zi(t)). Using this formulation the control dynamics equations to 
solve are: Eqs. (f|), (EDD) dUD with the boundary conditions Eq. fll5l) . This approach will 
be called method II. 

We can reformulate another cost functional based on replacing Eq. (0) with 

d(Zi(t)) 



dt 



/i((z(t)),e(t)) i = l,...,2N (20) 



This equation signifies that the time dependence of the expectation values (the quantum 



motion behavior) resembles one classical trajectory in the 2N dimensional quantum phase 
space. The cost functional may be now rewritten as 



T 27V 



J[(*),e] = J[(*),e] - I dtJ2\(t) 

Jo ~1 



(21) 

i=i L Ui J 
The variations of J[(z), e] with respect to (zi(t)), e(t) and X(t) produce the control equations. 

The variations with respect to Aj(i) give Eq. (|20|) . The variations with respect to (zi(t)) 

gives the same boundary condition as Eq. fll5D . The equations of motion for the Lagrange 

multipliers are identical to Eq. ([19]). However, the gradient of the cost functional with 

respect to electric field is different, 

6Jl(z),e] _ diMt)) . ™.,.. dfmt))At)) 

These control of equations are coupled to the Schrodinger equation. This approach will be 
referred to as method III below. To solve this system of non-linear differential equations we 
need to know the expectation value for each observable in phase space. The quantum control 
feedback equations (|I5D, ( |T9| ) and Eq. (^) are identical to classical feedback dynamics with 
a single trajectory, the average trajectory; M a comparison between these quantum and 



classical control dynamics equations will be reported elsewhere. |21| 

To summarize, the conventional approach in II. A and methods I, II and III in II. B should 
be viewed as providing alternate routes to equivalent control field designs consistent with 
the proper quantum dynamics of the molecule and the physical objectives. 

III. SELECTIVE CONTROL OF A UNIMOLECULAR REACTION 

The conventional quantum optimal control dynamics equations Eqs. (|j), (H), © and (§) 
in section II. A has been applied to a variety of problems including selective bond breaking 
through infrared excitation. In order to compare with the approaches described in Section 
II. B, we treat the same model expressed previously. The selective dissociation of one bond 
in a linear triatomic molecule is studied. The molecule is modeled as a pair of kineticly 



coupled Morse oscillators. |23[ The molecular Hamiltonian in internal coordinates is 



H mol = P- + p - 4r P ^ + V ( R i) + V(R 2 )mol (23) 

l\X\ Z/I2 M C 

where 

V(R t ) = D l [l-exp(-a l R t )] 2 , (24) 

R 1 , R 2 are the displacement operators of the bonds from their equilibrium positions, P 1 , P 2 
are the conjugate momentum operators corresponding to R± and P 2 , the reduced masses 
(amu) are ^ = ^ + ^ ^ = ^ + M c = 12amw, and 



P>i = 0.1976aw, D 2 = 0.13705aw, 
a\ = 0.9217aw and a 2 = 1. 02725cm. 
The dipole function of the molecule has the form 

R 2 ) = /i e [(i?! + R ) exp {-PRi) - (R 2 + R ) exp {-(3R 2 )} (25) 

where the parameters are \i e = 0.3934aw, R Q = 2au and (3 = law. The polarization of the 
external field e(t) is along the axis of the molecule. 

The model treats the selective dissociation of bond 1 at a minimum disturbance of the 
non selected bond 2 along with a minimum fluence of the electric field. Following the early 
studies (Ref. [2]) we choose the objective function to be 

*[<3(T)>] = \[(Ri(T)) - l] 2 P n + i(Pi(T)) 2 M-(P 1 (T))]P /3 (26) 

where 7 is the "target stretch", Pf 1: Pj 3 are positive constant weights and h(x) is the 
Heaviside function 

ri, ifx>o 

h(x) = (27) 
[0, if x < 0. 

The constraint function £i(z(t)) is chosen as 

4(z(t)) = i[(i? 2 (t)> - <P 2 (0)>] 2 ^ 2 + \[{P 2 (t)) - <P 2 (0)>] 2 ^ 4 (28) 
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where W 2 ,W4 are positive weights. This function biases the dynamics such that bond 2 
remains minimally excited and (i?2(t)) does not move far away from its initial value. The 
last function £2(e(t)), is chosen to minimize the intensity of the laser pulse 
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(29) 



where u e is a positive weight. 

Utilizing the above criteria for the cost functional Eq. (|TT|), where z = [zi,Z2,z 3 , z<±] repre- 
sent the coordinate and momentum operators, we obtain the corresponding set of equations 
(|18|) for the Lagrange multipliers of approach I 
dXx(t) 



dt 



As(f) 



(30a) 



d\ 2 (t) 
dt 



m)\v\z 2 )\m)+^m)\^'{z 2 )m)W) 

-W 2 [(z 2 (t)) - (z 2 (0))] 



(30b) 



dXi(t) 
dt 



^ = -M!) + J_ Aa(t) 

dt fjLi M 



■^- + ^Xi(t)-W 4 [(z 4 (t))-(z,(0))} 



where the prime denotes the derivative with respect to Zi, e.g. V" (zi 
boundary conditions Eq. fll5]) are 



A 1 (T)=P /1 [(2 1 (T))- 7 ] 
A 2 (T)= 



A 3 (T)= P f3 (z 3 (T))h(-(z 3 (T))) + ^(z 3 (T)y5(-(z 3 (T)}) 
A 4 (T)= 



(30c) 



(30d) 



qpl. The 



(31a) 
(31b) 
(31c) 
(31d) 



and the gradient of the cost functional J[(z), e] with respect of the field in Eq. ( |17D is 
5J[( Z ),e] 



5e(t) 



tv e 6(t)+Lie[Ht)mw(zi)\m) - \4(t)(m\Az2)\m)} 



(32) 
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The minimization condition SJ g^^ = gives 

ei {t) = -^Mt){mw(*i)\m) - wwwfamt))] (33) 

The electric field depends on the two Lagrange multipliers Xs(t) and X^it) and the average 
value of the dipole function derivative. Equation (|33|) links the Lagrange multipliers with 
the Schrodinger equation of motion. 

Also for method II we treat the Lagrange multipliers in the same manner as in Eq. ( |T§D . 
These equations are 

dXi(t) 



dt 



v (( Zl (t))) - w"(( Zl (t)Mt) 



(34a) 



dX 2 (t) 
dt 



A 4 (t) 



V"{(z 2 {t)))+^"{(z 2 {t)))e{t) 



\{z 2 {t)) - (z 2 (0))]W 2 (34b) 



dX 3 (t) XAt) 1 , . . . , 

-5T = "TT + M h(t) (34c) 

^ = ~ + Tf Xl(t) ~ l<Z4(i)> " {z ' mWi (34d) 

The control equations Eq. fl34 ) have the same boundary conditions Eq. (|31~1) and the same 



gradient of the cost functional with respect to electric field in Eq. ( J32|) . In this approach II 
the Lagrange multipliers will differ from those of the previous approach I and therefore the 
minimum solutions €j(t) and e/j(i) are expected to be different. 

The approach III has the same Lagrange multiplier equations as in the approach II in 
(j34|) and the gradient of the cost functional Eq. (|21|) with respect to the electric field is 

6J[(z),e) 



5e(t) 



u e e(t) + // e [A 3 (*y - M*)/A<*2(*)»] (35) 



The minimization condition SJ ^^ = gives 



ein{t) = -^[X 3 (t)^(( Zl (t))) - X 4 (t)»'((z 2 (t)))} (36) 
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Note that these equations are identical to those reported previously |T| for molecular clas- 
sical optimal control evaluated over the average trajectory. These set of quantum control 
equations Eq. (|3~3|), fl35|) with the boundary conditions Eq. fl3~T|) can be treated consistently 
with the Schrodinger equation. Solving the Schrodinger equation directly, the trajectory of 
the expectation values (zi(t)) are obtained and are then utilized to solve the coupled set of 
equations Eq. (0), (|35|)- This is the basis for the self consistency of this approach. 

In the following section we explain the computational method used to solve the equations 
leading to e/(t),e//(t), e HI {t). 



IV. COMPUTATIONAL METHOD 

To explain the computational methodology first consider the approach I. The set of control 
dynamics equations to solve are Eqs. (|3|), (|30|), ([32]) with the boundary conditions Eq. (PH). 
An iterative scheme is adopted with the aim of finding ej(t) at a minimum cost while meeting 
the desired objective. The control algorithm used in this paper is the following: 

a) Make an initial guess for the field ei{t). 

b) Integrate Schrodinger's equation Eq. (|3]) and calculate the quantum average trajectory 

and other expectation values necessary to integrate Eq. (p0|). 



c) Calculate the probability of dissociation for each of the channels and evaluate the cost 
functional. 



d) Integrate Eq. (|30D for the Lagrange multipliers backward in time. 

e) Calculate of the gradient of the cost functional with respect to the electric field. 

f) Update a new field from the last step. 

The steps (b) to (f) are repeated until a converged solution is obtained. 
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Similar algorithms can be constructed for the two other approaches II and III. In order to 



improve the form of e(t) on each iteration we use the conjugate gradient algorithm. p0| 



The Schrodinger equation was integrated using the split operator method. flT5l , |I~4|| This 
requires the propagator 

rt+St 

it 



U(t + St, t)« exp(— / [V(R) + H int (R, e{t'))]dt'/h) exp(-iK5t/h) 
2 Jt 
_2 rt+St 

exp(— J [V(R) + H int (R, e(t'))}dt' /%) (37) 



which evolves the wave function St in time. The number of grid points chosen for the two 
coordinates [Ri^R^) was 256 x 128. The range of values for R\ and R2 was [-l,12]aw and 
[-l,5]cra, respectively. The time step St = 16. law was used to propagate the wave function 
over a time interval of O.lps, and for times larger than O.lps, St = 0.5cm was selected. 

An optical potential was introduced to absorb flux in the boundary region. The new 
Hamiltonian is 

H = H mo i + H int — iVop t (R) (38) 

where the optical potential is chosen to be 

(V %=z, if a <Ri<b; 

= " (39) 

[ 0, otherwise. 

The parameters for the first bond are V a \ = 2a. u., a = 11.5a.u,b = 12.0a.-u and those for 
the second bond are V Q 2 = la.u, a = 4.5a.u, b = 5.0a.u. 

With this scheme, we can define the dissociation probability in various ways. The most 
straightforward approach is simply to compute the probability as 

Pi(t) = f dt' t dR d Ji(R d , t') (40) 

JO Jen 

where the J~i(Rd,t) is the flux defined as 

Ji(R d ,t) = -Im(tf>(t)\h-Vi\^{t)) (41) 

The spatial integral in Eq. QlTj) is over the flux dividing line between two points in the 
surface a« and hi, and n is the local unit vector normal to the dividing line. The gradient 
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in Eq. (f41~|) is with respect to Ri, R2 and the index i refers to a specific channel in the 
photodissociation process. We can get four different channels in the process of dissociation, 

' A + BC dissociate bond 1 

AB + C dissociate bond 2 



ABC 



(42) 

A + B + C full break up 



- ABC no dissociation 

where each process has a distint probability of occurrence. In order to calculate the proba- 
bility for each channel, we choose the following values: Pi = A + BC for the flux through 
the line [6,-1] to [6,4]; P 2 = AB + C for [-1,4] to (6,4); P 3 = A + B + C for [7,4] to [6,5] (all 
of these values are in a.u.). The probability for P4 = ABC was calculated as 

P<(t) = 1-1; Pi(f). (43) 

i=i 

V. NUMERICAL RESULTS 

In Fig.l we show the optimal laser pulse and its corresponding power spectrum, utilizing 
the standard quantum control equations Eqs. (|3|), and (§) to (^j) of section II. A. The 
optimal laser pulse is identical to that in Ref. [2]. The power spectrum of the field has only 
one dominant peak at 1334cm -1 , which is higher than the two fundamental frequencies 
u° = 1098cm -1 and u)^ = 923cm" 1 of both Morse oscillators. The second peak at 667cm _1 
plays a less important role. 

In Fig. 2 we show the optimal pulses obtained from the three different algorithms I, II 
and III proposed where the guessed input field was zero. The optimal pulse I is produced 
from the solution of the control dynamics equations (|3|), fl30| ) and ( |32|) with the boundary 
condition Eq. (pT|). The optimal solution II is obtained with the Lagrange multiplier given 



by Eq. ([34]). The result III in Fig. 2 is based on the control dynamics equations calculated 
over one quantum trajectory utilizing Eqs. (||), fl3"4|), (ft5|) and Eq. (|3"T|). In all three cases 
the optimal solutions for the electric field selectively dissociate bond 1. All three fields ej, 
ejj and em are strikingly similar and also closely like that in Fig. 1. The power spectra of 
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ej, en and em are also similar to that in Fig. 1, except that the small low frequency peak 
near 700cm" 1 is missing. 

Figure 3 shows the temporal evolution of the quantum expectation values for the approx- 
imate bond energy, bond lengths, momenta and total molecular energy in the presence of the 
optimal field of Fig. 1. The analogous results for the field ej of Fig. 2 is shown in Fig. 4 (the 
results of applying en and em are almost the same as ei). The desired dissociation event 
ABC^ A + BC is clearly evident. The approximate energy for each bond (this energy is 
approximate because the kinetic coupling between the two bonds is not taken into account) 
shows high excitation in bond 1 (the dissociation energy is 0.197a. it) much in contrast to 
bond 2 (the dissociation energy is 0.137a.w). The time evolution of the expectation values 
for the positions and the momenta for the two bonds also indicate selective dissociation 
of bond 1. The expectation values of the interaction energy and the total energy of the 
molecule as a function of time are shown in Fig. 5 for the conventional field in Fig. 1 and 
the approach I (the results of II and III are essentially the same as the latter case). We 
observe that the interaction energy has its maximum value around the time t = 0.055ps in 
all the cases. 

Table 1 shows the probability for the four reaction channels using the four optimal pulses 
of Figs. 1 and 2. More than of 50% of the dissociation occurs for the desired channel ABC 
— > A+BC with less than 1% for the channel ABC — > C+AB in all cases. The low values 
for the dissociation probabilities for the reactions ABC — > C+AB and ABC — > A+B+C 
demonstrates the high degree of control achieved for the new approaches to designing controls 
over quantum motion. The number of iterations to achieve the same level of convergence with 
the various methods is essentially the same as shown in table 2. However, the computational 
saving of methods I, II or III over that of conventional quantum control is nearly a factor of 
two as the overhead is very small for solving the Lagrange multiplier equations with the new 
methods. The potential for further computational saving will be discussed in the conclusion 
section. 



16 



VI. CONCLUSIONS 



In this work we have presented a new general formalism for the control of quantum ob- 
servables. Illustrations were carried out to compare a family of related methods. The 
conventional reference approach standardly utilizes a full Lagrange multiplier state vector 
introduced to preserve the Schrodinger equation. The Lagrange multipliers in the new 
approach are scalars introduced to preserve the average dynamics. A family of related ap- 
proaches I, II, III was presented, where each member corresponded to a distinct treatment 
of the Lagrange multiplier equations. 

The four algorithms presented in this work were applied to a triatomic molecule model 
studied previously. The numerical results demonstrate the selective dissociation of the 
stronger bond 1, while bond 2 remains only weakly excited. The approximate bond en- 
ergy and the temporal evolution of the bond length and momentum expectation values 
demonstrate selective dissociation for each of the optimal fields. We also compared the ex- 
pectation values of the total molecular energy with the interaction energy as a measure of 
the efficiency of the control process. From these observations we deduce that the three new 
algorithms I, II, III proposed here give excellent and very similar results. The dissociation 
probability of the desired bond 1 was more than 50%, while for bond 2 the dissociation was 
less than 1%. 

With these excellent results in evidence, a central question is why such a serious alteration 
of the Lagrange multiplier equations is successful. Three factors are operative here: 

(1) The proper quantum dynamics of the molecule through \ip(t)) is fully retained. 

(2) The cost functional retains the goal of achieving the original target objective. 

(3) There are typically multiple (if not an infinite number) of possible control solutions in 

any physically well posed molecular control problem. 

The role of the Lagrange multipliers is to provide feedback and guide the dynamics to an 
acceptable solution. Considering the three points above, the alternative formulations we 
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introduced for \\(t)) just serve to take us from one acceptable solution to another. The 
striking similarity of the reference pulse and those of methods I, II and III also suggests 
that the field minimization space of the cost functional is not locally distorted to a signif- 
icant degree. The ultimate saving with the new methods resides in their simplification of 
the traditional strongly coupled two point boundary value problem for \4>(t)) and |A(t)). 
The alternate methods for |A(t)) greatly reduce the complexity of this task. The different 
dynamics for X(t) in methods I, II and III give almost the same fields (actually, fields II 
and III are numerically the same) which in turn are very similar to that of the conventional 
approach. This strong similarity may not always occur in other problems, but its presence 
here clearly indicates the wide latitude in treating the feedback process in control field de- 
sign. Capitalizing on this flexibility, by the approaches suggested here or other related ones, 
could greatly simplify the molecular control design process. 

In the present paper the new Lagrange variable acted to preserve the dynamics for the 
molecular bond length and momentum expectation values. The same logic could also be 
applied to preserve additional or distinctly different observable expectation values besides 
(zi(t)). Analogous problems of this type also arise in more traditional engineering applica- 
tions of control theory, and it would be interesting to apply these reduced Lagrange multiplier 
concepts in such cases. 

The ultimate significance of the findings in this paper is suggested by considering the 
work in the context of the analogous classical study |J and recent tracking control stud- 
16Hl9f Tracking is relevant here, as it operationally replaces the feedback role of the 
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Lagrange multiplier by an expression for the control field explicitly in terms of the system 
wave function. In essence, the present paper introduces what may appear to be serious 
approximations for the feedback Lagrange multipliers; but, in fact, the resultant control is 
well-achieved, and in some cases, the field is strikingly similar to that obtained by the full 
traditional feedback approach. Similarly, the operations of tracking would appear to create 
a drastic modification of traditional feedback, yet tracking encompasses traditional optimal 
control solutions, as well as others. All of this work points to the observation that there is 



considerable freedom for introducing approximations and direct physical guidance into the 
feedback aspects of quantum control. This feedback can be compactly expressed as 

Am) 



ih- 



H mol + ne(t, (0),\iP(t))) 



\m) (44) 



dt 

where the field is shown to possibly depend on time explicitly the expectation value (O) of 
target observable operator O and the state \ip{t)) of the system. Tracking control has this 
form of feedback, and a formal solution to the Lagrange multiplier equations also leads to 
a similar form of feedback. In the simplest case, Eq. fl44j) only needs to be integrated once 
to achieve a design. The ultimate savings from these overall simplified approaches could 
be a factor of two, or up to many times that magnitude, when iteration is eliminated as 
in tracking and other direct feedback approaches (i.e. for tracking the factor of savings is 
~ 2N where N is the nominal number of optimal control iterations and tipically N » 10). 
The insight gained from the present body of work suggests that focusing on the physical and 
numerical content of the feedback should be a very fruitful direction in molecular control 
theory. 

This new approach may be conveniently combined with methods that compute the po- 
tential surface as the dynamics proceeds. [P9ppE5|] The Lagrange multiplier logic has now 



been applied to classical mechanics [|l],|4|-|26[ and quantum mechanics in the present paper. 
The same control concepts may also be applied for semiclassical wave packet propagation 
HJ-46 and mixed quantum/classical molecular dynamics. |35|-p8 
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FIGURES 

FIG. 1. The optimal electric field and power spectrum for selective dissociation using the 
traditional approach of Eqs. (||), @, (f|) and Eq. @. The system is in the ground state at 
t = 0. The parameters in the control dynamics equations are 7 = 5, Pfi = 400-Diaf, Pf% = 222, 
W2 = 2D2Q.2, W4 = -j^. The pulse duration time was O.lps. These parameters also apply to the 
following figures. 

FIG. 2. Optimal electric fields: Method I using the Eqs. @ to @); Method II using the 
Lagrange multiplier in Eq. (j34|); Method III using the Eqs. ( [34] ) and (|35| ) with the same boundary 
conditions. The results of methods II and III are numerically identical. 

FIG. 3. Time dependence of the quantum expectation values of the approximate bond energy, 
bond length and momentum for each bond for the optimal pulse of the Fig. 1. 

FIG. 4. The same as Fig. 3 for the optimal pulse I of Fig. 2. The results from the pulse II and 
III are similar to those of pulse I. 

FIG. 5. Time dependence of the quantum expectation value of the total molecular Hamiltonian 
(Ex) and the interaction term (Ei nt ). a) for the reference pulse in Fig. 1; b) for the pulse I of Fig. 
2. The results from the pulse II and III are very similar to that of pulse I. 
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TABLES 

TABLE I. Convergent value of the dissociation probability Pj for the four reaction channels. 
The four channels are: Pi for ABC^ A + BC; P 2 for ABC^ C + AB; P 3 for ABC^ B + A + C 
and P4 means no dissociation. CM means the conventional method (the pulse of the Fig. 1). For 
the approaches II and III we obtain almost the same dissociation probability as that in approach 
I. 



Method Pi P 2 P 3 P 4 

CM 0.560 1.8xl0 -3 lxl0~ 3 0.437 

I 0.628 7.3xl0~ 3 2.7xl0~ 3 0.362 



TABLE II. Number of iterations for the four different methods in order to find the optimal 
pulse field in Fig. 1 and Fig. 2. 



Method Number of Iteractions 

CM 17 

I 18 

II 19 

III 15 
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Figure 1 
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Figure 4 




Figure 5 




